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ABSTRACT 



The Quintuplet, one of three massive stellar clusters in the Galactic center (GC), is located about 30 pc in projection from Sagittarius 
A*. We aim at the construction of the Hertzsprung-Russell diagram (HRD) of the cluster to study its evolution and to constrain its 
star-foiTnation history. For this purpose we use the most complete spectral catalog of the Quintuplet stars. Based on the A"-band spectra 
we determine stellar temperatures and luminosities for all stars in the catalog under the assumption of a uniform reddening towards 
the cluster. We find two groups in the resulting HRD: early-type OB stars and late-type KM stars, well separated from each other. By 
comparison with Geneva stellar evolution models we derive initial masses exceeding 8 Mq for the OB stars. In the HRD these stars 
are located along an isochrone corresponding to an age of about 4 Myr. This confirms previous considerations, where a similar age 
estimate was based on the presence of evolved Wolf-Ray et stars in the cluster. We derive number ratios for the various spectral subtype 
groups (e. g. A'wr/A'o, A'wc/A'wn) and compare them with predictions of population synthesis models. We find that an instantaneous 
burst of star formation at about 3.3 to 3.6 Myr ago is the most likely scenario to form the Quintuplet cluster. Furthermore, we apply 
a mass-luminosity relation to construct the initial mass function (IMF) of the cluster. We find indications for a slightly top-heavy 
IMF. The late-type stars in the LHO catalog are red giant branch (RGB) stars or red supergiants (RSGs) according to their spectral 
signatures. Under the assumption that they are located at about the distance of the Galactic center we can derive their luminosities. 
The comparison with stellar evolution models reveals that the initial masses of these stars are lower than 15 Mq implying that they 
needed about 15 Myr (RSG) or even more than 30 Myr (RGB) to evolve into their present stage. It might be suspected that these 
late-type stars do not physically belong to the Quintuplet cluster. Indeed, most of them disqualify as cluster members because their 
radial velocities differ too much from the cluster average. Nevertheless, five of the brightest RGB/RSG stars from the LHO catalog 
share the mean radial velocity of the Quintuplet, and thus remain highly suspect for being gravitationally bound members. If so, this 
would challenge the cluster formation and evolution scenario. 
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1. Introduction 

Stellar clusters are unique objects to study stellar populations 
and their evolution, as they are supposed to represent a single co- 
eval stellar population. The Quintuplet cluster is one of the mas- 
sive young stellar clusters in the Galactic center region. Located 
at about 30 pc projected distance from the Galactic Center, its 
age has been estimated to be about 4 Mvr (iFiger et al.l 1 1 999al) 
from the evolved massive stars in their Wolf-Rayet (WR) phase, 
i.e. stars displaying CNO-processed (WN stars) or helium burn- 
ing products (WC stars) in their spectra. T he WN stars were an- 
alyzed in detail bv lLiermann et al.l(l2010h who found from com- 
parison with stellar evolution models ages of about 3 Myr that 
agree quite well with the previously assumed cluster age. 

Despite the young cluster age, iGlass et aTl (Il990h reported 
the presence of one red supergiant (RSG) in the cluster A 
few other Galactic clusters are known to show the coexistence 
of evolved stars like WR stars and RSGs, e. g. Westerlund 1 
(ICrowther et al.ll2006l) . but it seems quite puzzlin g for a young 
cluster like the Quintuplet. iLiermann et al.l(l2009l) presented the 
so far most complete spectral catalog of the Quintuplet stars. 
From this catalog the Heitzsprung-Russell diagram (HRD) can 
be constructed, and the results can be compared to stellar evolu- 
tion models and isochrones. This will allow to put more stringent 



constrains on the age of the Quintuplet cluster The comparison 
of number ratios with population synthesis models will help to 
confine the cluster age further and allows a statement about the 
formation history to be made. Especially, the apparent existence 
of RSGs in the cluster needs to be addressed. 

The paper is organized as follows. In Sect.|2]we introduce 
the spectra and the sample of stars, Sect.|3]explains the methods 
to deiive a cluster HRD. Finally, the results are discussed in the 
context of stellar and cluster evolution in Sect.|4] 



2. Observations 

The analysis in this paper is based on observations of the 
Quintuplet cluster with the ESO Very Large Telescope's 
integral-field spectrograph SINFONl, that cover /T-band spectra 
with a resolution of about R x 4000, complete to a photometric 
magnitud e of about K^ = 13 mag. The data were published in a 
catalog bv lLiermann et al.l (120091 hereafter LHO catalog), listing 
85 early-type and 62 late-type stars with their spectral classifica- 
tion and individual radial velocity (RV). A field of view of about 
40 "x 40 "is covered in total. Please refer to the LHO catalog for 
more detailed information about observational strategy and data 
reduction. 



A. Liermann et al.: The Quintuplet Cluster. III. 



50 

40 

z: 
3 30 




la 



3a 



40 80 120 

radial velocity [km/s] 

Fig. 1. Histogram of the distribution of /JVs of all LHO stars 
(light gray). Overplotted is the distribution for the late-type KM 
stars (dark/red). Solid and dashed lines indicate the cluster RV 
of 113 km s"' and the 1 cr and 3 cr intervals, respectively. 
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Fig. 2. Same as Fig.lT] but overplotted with the RY distribution of 
early-type OB stars (dark/blue). Both distributions peak within 
the intervals assumed for cluster membership, albeit more pro- 
nounced for the early-type stars. 



3. Analysis 

3A. Sample stars 

The LHO catalog gives a mean cluster radial velocity of RV - 
113 + 37kms"' and a first approximation for the velocity dis- 
persion cr - 17kms ' within the Quintuplet cluster, under the 
assumption of a virialized cluster The cluster membership was 
rejected for stars with a radial velocity that differs from the clus- 
ter mean by more than 3 cr (+ 51 kms"'). However, field stars 
in the vicinity of the Quintuplet may participate in a similar 
Galactocentric rotation and thus show a similar RV, albeit not 
being bound to the cluster Moreover, RVs can be affected by the 
orbital motion for stars in close binary systems. In these cases 
the periodic variation of the RVs can only be detected from time 
series of spectra, which are not (yet) available for the Quintuplet 
stars. 

The number distributions of the radial velocities of the late- 
type and early-type stars are shown in Figs.[T]and |2] While we 
find a very pronounced peak in the distribution for the early- 
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Fig. 3. Distribution of Ks magnitudes over radial velocity for all 
stars in the LHO catalog (circles for early-type stars, triangles 
for late types). The completeness limit at Ks - 13 mag is indi- 
cated by the horizontal line, while the vertical lines refer to the 
different RV intervals as in Fig.[T] 



type stars in the RV range assumed to indicate cluster member- 
ship, the late-type KM stars seem to be distributed more homo- 
geneously with a somewhat lower peak within the 3cr interval 
around the mean cluster RV. To estimate the number of stars, 
late-type stars in particular, along a line-of-sight towards the 
Galactic center, we would need to observe control fields around 
the Quintuplet. The lack of those requires a more careful inter- 
pretation of Fig.[T] from the rather flat distribution, we could as- 
sume a median number of five stars per velocity bin to possibly 
be field stars. 

Additionally, the question of foreground stars has to be re- 
garded under the aspect that the high visual extinction towards 
the Galactic center is overcome in the near-IR range. This in re- 
turn makes it more difficult to distinguish between foreground 
and background objects relative to the cluster As can be seen 
from Fig. |3] there is no clear distinction in magnitude between 
bright foreground or faint background objects, but a wide distri- 
bution of the TiT-band magnitudes over the RV range. The photo- 
metric completeness of the LHO catalog at Ks - 13 mag hardly 
alters the distribution. 

We set the following criteria to assess the possible clus- 
ter membership: stars within the 1 cr interval of the cluster RV 
and brighter than the photometric completeness limit, i. e. Ks < 
13 mag, will be considered as very likely cluster members. We 
will refer to this sample of stars as the corrected sample in the 
following in comparison to the complete LHO sample. Q 

3.2. Stellar luminosities and the HRD 

The HRD is a useful tool to characterize a stellar cluster; to es- 
tablish it we need to know the effective temperatures and stellar 
luminosities for the Quintuplet stars. 

The effective temperature is defined via the Stefan- 
Boltzmann law, and thus follow from the stellar luminosity and 
radius. This reference radius is a matter of definition. Usually 
one refers the stellar radius to a Rosseland optical depth of 2/3. 



' The criteria are slightly stricter than the previous 3 cr cluster RV, 
also applied in the LHO catalog, but are meant to compensate to some 
degree the lack of time series of spectra to exclude binary effects and 
the lack of control fields to estimate the number of field stars properly. 
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In stars with a well defined photosphere the precise definition 
does not matter In WR atmospheres, however, t - 2/3 is often 
reached already in moving layers, i. e. in the wind. Therefore it 
became standard to define the "stellar radius" /?» of WR stars at 
a Rosseland optical depth of 20. In most models, this point lies 
at subsonic velocities, i.e. in nearly hydrostatic layers. The effec- 
tive temperature related to /?, is termed the "stellar temperature" 
r», and this value is quot ed for the WN stars as the result from 
the spectral analyses from lLiermann et all (1201 Ob . 

The motivation for using the radius of the hydrostatic part 
is the hope that /?» can be identified with the stellar radius from 
stellar evolution models. However, our empirical R, is in fact 
based on an inward extrapolation of the stellar wind's velocity 
law into optically thick, i.e. un-observable layers. In the case 
of our Quintuplet WN stars, this is not critical because R, is 
only slightly smaller than the radius where t - 2/3 - only 
in very thick WR winds, the difference becomes significant. 
Nevertheless, there is evidence that most WR stars have very ex- 
tended subphotospheric layers which make /?* much bigger (and 
thus r* much lower) than current models of the stellar structure 
predict. 

For the WC stars in the Quintuplet the analysis is in progress . 
All o f them are classified as WC8-9 and their m ajority produces 
dust dvan der Huchtll2006l iLiermann et al.ll2009l) . This dust pro- 
duction is considered to arise from colliding winds in high- 
mass binaries, as has b een confirmed alread y for some of the 
Quintuplet WC stars by iTuthill et all (l2006i "pinwheel stars"). 
Therefore, the WC stars will be excluded from the present con- 
siderations of the cluster HRD, as they need to be analyzed and 
discussed with respect to the special scenarios of binary evolu- 
tion. 

For the late-type KM stars, we determine the effective tem- 
peratures (Teff) according to their spectral classification in the 
LHO catalog. Those were based on the equivalent widths mea- 
sured from the first overtone CO absorption band at 2.3 /vm (see 
iGonzalez-Fernandez et al. 2008) . The Teff are then applied to ob- 
tain the bolometric correction for the /T-band magnitude, BCa:, 
adopting the relation 

BCk = 5.574 - 0.7589 (r^ff/lOOOK) , (1) 

derived bv lLevesque et al.l (l2005h for red supergiants. We apply 
the same relation for our red giant branch stars (RGBs) as well, 
thus neglecting effects of the luminosity class on the spectral 
energy distribution. 

The effective temperatures for the early-type OB type stars 
were determined from their spectral (sub-)type s accordin g to the 
spectral-type-temperature calibration of Mar tins et al.l (I2005L 
Tables 2,3,5, and 6) and applied for their bolometric corrections 



7 - 



BCy = 27.58 - 6.80 log (T^s/K) . 



(2) 



These BCs are valid only for the visual spectral range as in- 
dicated by the index and a correcting term for the /T-band 
BCs is necessary. From a previous model analysis of O stars 
JOskinova et al.ll2006.) we can estimate an average My - Mk — 
(-0.98 + 0.04) mag for their sample stars. That sample contained 
spectral types 03 to 07 with luminosity classes I to V, there- 
fore the term is considered to be robust and applicable for the 
correction BCa: - BCv - {My - Mk)- 

Absolute stellar magnitudes Mk were derived from the K^ 
given in the LHO catalog with the mean cluster extinction 
Ak - 3.1 + .5 mag derived from the analyses of the WN stars 
(ILiermann et al. 201 0) together with the corresponding BCs. For 
the cluster distance we adopt the distance to the GC of 8kpc 
(lReidll993h . 
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Fig. 4. Hertzsprung-Russell diagram of the Quintuplet. Circles 
(blue) represent the early-type OB stars, triangles (red) the late- 
type KM stars, filled symbols refer to stars of the corrected sam- 
ple (see text). The ZAMS and stellar evolution tracks with rota- 
tion (solid lines) and without rotation (dashed lines) for different 
initial masses are from Mevnet & Maeder (2003). 



The absolute magnitudes are transformed to stellar luminosi- 



ties 



log — = -0.4 (M;f - BC^ - Mbouo) , (3) 

with Mboi,0 = 4.74 mag for the bolometric luminosity of the Sun. 
The resulting HRD of the Quintuplet cluster is shown in Fig.|4] 

The errors on the temperature and lu minosity of the WN stars 
are taken from lLiermann et all (1201 Ob . 

For the late-t ype KM stars, the e rror on T^^ is of the order 
of 200 K (Gonzale z-Fernandez et al.l f2008). For those stars the 
error bars vanish within the stars' symbol sizes. 

For the early-type OB stars the error on Teff reflects the range 
in the spectral classification in the LHO catalog as an upper and 
lower temperature according to the spectral-type-temperature 
calibration. This error is usually quite large and could be min- 
imized by tailored modeling of the individual stars, which is out 
of scope for the work presented here. 

For both late- and early-type stars we estimate the error on 
the luminosity from the combined errors of the A'-band magni- 
tude, extinction, and bolometric corrections. The error is dom- 
inated by the error of the mean cluster extinction determined 
from the WN star analysis (0.5 mag). The error on the /T-band 
mag nitudes from the LHO catalog is given with 0.2 mag. For the 
BCs Levesque et al. (2005) list an error of 0.01 mag for late-type 
stars and Martins et al. (2005) give 0.05 mag for early-type stars. 
We had to correct the BCs for the early-type stars from the vi- 
sual to the /T-band with a correcting term for with we estimate 
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Fig. 5. Same as Fig. |4] but with theoretical isochrones for dif- 
ferent cluster ages. The isochron es, shown with alternating line 
styles, were constructed by Leieune & Schaerer (2001 ) . The ma- 
jority of the population of OB stars in the Quintuplet follows the 
4Myr isochrone, while most KM stars would need more than 
SOMyrto evolve 



an error of 0.04 mag (see above). In total, the combined error on 
the luminosity amounts to about 0.3 dex for both the late- and 
early- type stars. 



4. Discussion 

4.1. Stellar evolution and cluster age 

As can be clearly seen from Fig.|4] two distinct groups are found 
in the HRD. The OB stars are located more or less along a main 
sequence, while the KM stars line up in the low-temperature 
regime with Teff < 5000 K over an extended luminosity range. 
It seems that a number of OB stars are found to the left of the 
zero-age main sequence (ZAMS). On one hand, this can be at- 
tributed to the rather large error bars in temperature, i.e. uncer- 
tainty of the star's spectral type, and luminosity. However, it 
might also reflect the fact that the assumed homogeneous red- 
dening towards the cluster is not correct. And with the presence 
of dusty WC stars an additional intrinsic reddening might have 
to be considered for the diff'erent regions of the cluster 

Stellar evolution tracks including the effects of rotation (solid 
lines) and without r otation (dashed Unes) for initial masses of 9, 
15, 25 and 60 Mq dMevnet & Maedeill2003h show that the OB 
stars are massive stars in the classical definition, i. e. Minit > 
8 Mq. We omit the tracks for higher masses as they overlap with 
the track for 60 Mq in case of the 85 Mq track or don't extend 
into the cooler temperature range in case of the 120 Mq track. 
The difference in stellar evolution between models with rotation 
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Fig. 6. Same as Fig . |5] but with theoretical isochrones based on 
iGirardi et al.l (12002). Again, the majority of OB stars falls on the 
4Myr isochrone, while most KM stars evolve on timescales with 
more than 30 Myr 



and without is stronger from 25 M© upwards, i. e. the high-mass 
stars potentially evolving to WR stars, and is less important for 
the majority of stars in the current sample. 

Concerning their position in the HRD below the 9 Mq track 
and off the main sequence, the late-type KM stars would have 
to be considered evolved low- and intermediate stars. This is 
confirmed by the presence of ^^CO absorption in their spectra 
which led to the classification as (super-) giants (LHO catalog). 
However, a few late-type stars extend to the regime of massive 
stars touching the 9 to 15 Mq tracks. Their potential status as red 
supergiant (RSG) will be discussed separately (see below). 

In Fig.[5] the HRD with isochrones compiled by 
iLeieune & Schaeren (12001, based on Geneva stellar evolu- 
tion models) is shown for different cluster ages . For co mparison 
we show the isochrones from |Girardi e t al.l (l2002[ "Padova 
isochrones" in the following) in Fig.|6] In both cases, the as- 
sumed cluster main sequence of OB stars corresponds well with 
the 4 Myr isochrone with the 3 and 5 Myr isochrone forming an 
"age envelope". 

In case of the WN stars, two stars (WR 102d and WR 102i) 
seem to follow the Geneva 2 Myr isochrone, while the 2 Myr 
isochrone from Padova is below any WN star position (Fig.|6). 
In addition, both of these isochrones don't extend far enough into 
the cooler region of the HRD to cover all WN stars. Interestingly, 
the Padova isochrones for 3, 4 and 5 Myr loop backwards in the 
hotter region of the HRD, and especially the 3 Myr isochrone 
corresponds very well to the position of all the evolved WN 
stars. A detailed comparison of the W N stellar parameters with 
the Geneva stellar evolution models dMevnet & Maeden l2003h 
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Table 1. Number ratios for the Quintuplet stars in comparison to other Galactic center clusters. 
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Notes. Numbers of O stars are lo wer limits from iMartins et al.l (120081 Arches), iPaumard et al.l ( 120061 Central Parsec), iNegueruela et al.l d201j . 
Wd 1); WR numbers are based on Ivan der Huchj (l2006h unless otherwise indicated. WC stars showing dust in their spectra (WCd or WCLd 
classification) are considered binary candidates. Note that the number stars listed in the different columns for each clusters might slightly vary in 
dependence of the reference and the therein used subsample of stars (see text for details). 



* WR stars from lvan derHuchtJ (l2006h coiTected for the WR stars newly identified in lLiermann et al.l (l2009h : LHO 1 10 - WN9h, LHO 76 and 
LHO 79 both WC9 d; number of RSGs estimated from the corrected sample of late-type stars; five WCd stars resolved as pinwheel binaries by 
iTuthill p.t al. ^ ('2006), two further candidates based on WCd classification. 

** Numbers limited to the LHO field of view of about 40 "x 40 "coveri ng the center of the Quintuplet cluster. 

" Classification as WNL or WNE from Cot era et al.l (ll999) and Marti ns^t al.l l l2008h ; two WN bi naries from non-thermal, variable radio plus 
X-ray detection, further 5 candidates based on either variabl e radio flux or X-ray d etection dvan der Hucht 2006) . 

* Classification as \y NL or WNE from lPaumard et al.l(l2006h . lMartins et alJ ( l2008h : two Of/WN9 stars detected with radial velocity variations 
( IMartins et al.l2006h . 

'^ A^RSG includes four RSGs bv lMenge l & Tacconi-Gar manI (l2007l) and six yellow hypergiants bv lClark et al.l (l2005h : classification as WNL or 

WNE and WR binaries take n from'C rowther et al.l ( 120061) . 
" Numbers from lvan der Hucht (2001) . Table 29. A^vne^A^wnl ratio from lHamann et aP (l2006h . 



revealed ages of about 2.1 to 3.6 Myr for the Quintuplet WN 
stars (Lierma nn et al. 2010, 201 1 ) : an age of 3 Myr seems to be 
confirmed by the Padova isochrone. The presence of WC stars 
in the Quintuplet cluster speaks neither in favor of a much older 
nor a younger cluster age, since they may originate from binary 
evolution as explained above. 

The same argument holds for the luminous blue variables 
(LBVs) in the Quintuplet, the Pistol star and qF 362. From high- 
mass single star evolution the LBV phase can be estimated to 
start at abo ut 4 Myr, which was also used as age cr iterion for the 
Quintuplet (JFiger et al.ll999ailGeballe et al.l200"0l) . However, re- 
cent interferometric studies indicate that the Pistol star, just like 
the LBV prototype zyCar, might be a binary and thus evolved 
through the binary channel ("Martavan et al.l l201 ll) . 

A simil m' spread in ag e betw een O stars and WR stars was 
reported by iMartins et al.l (l2008l) for the Arches cluster They 
found ages of 2 to 3 Myr for the most luminous WN stars, while 
the O stars cover a range of 2 to 4 Myr They conclude that 
this might be due to the fact that the most massive stars have 
formed last in the cluster, to prevent their feedback on the on- 
going star formation. In the end they determine a cluster age of 
2.5 + 0.5 Myr Similar arguments can be applied for the the dif- 
ference in age of the OB and WN stars in the Quintuplet. Thus 
we conclude that a cluster age of 3.5 + 0.5 Myr is likely. 

Both sets of isochrones suggest that the KM stars, if they 
were cluster members, would need about 15 Myr to evolve to 
RSGs for Mini > 9 Mq and more than 30 Myr for Mini < 9 Mq to 
become red giants, respectively. The detected '■'CO gives spec- 
troscopi c evidence that these stars are evolved giants or super- 
giants (Liermann et al.l 12009.) . It appears that for the majority 
of these stars it might not be justified to assume that they are 
located at a distance of 8kpc. They can be foreground objects 
and/or physically unrelated to the Quintuplet cluster This con- 
clusion seems to be supported by the radial velocities of most of 
these stars which are off the limits of the corrected sample for 



the cluster However, we are left with a group of eight evolved 
KM stars that might be physical members of the Quintuplet (see 
below). 

A priori the coexistence of RSGs, main sequence OB stars, 
and WR stars is not expected under the assumption of a coeval 
evolution. However, iHunter et al.l ( l2000l) discuss the possibility 
of an extended star formation event lasting over a few million 
years, which could explain the simultaneous presence of WR 
stars and RSGs. We find an age range of about 4 Myr (OB stars) 
to 15 to 30 Myr (KM stars) for the two populations. Could such 
a spread in age be called an "extended" star formation event? 
Or is it more likely that one star formation event was responsi- 
ble for the WR and OB stars while the evolved KM stars might 
represent an older population that is unrelated to the cluster? 



4.2. Number ratios and the binary effect 

In Table[T]we list stellar number ratios for some prominent clus- 
ters and the Milky Way in total. These numbers are only ap- 
proximates. The ratios can suffer from systematic observational 
errors, incomplete samples used in the different references, or 
from unknown/unconfirmed binary stars in the samples. For ex- 
ample, among the general results of I van der Huchtl(l200ll) for the 
Galactic WR star population, they give a Nwc/Nwn ratio of 1.5 
and a binary frequency for the WR stars of 39%. These num- 
bers are based on a stellar subsample limited to WR stars within 
2.5 kpc; presumably a complete sample. On the other hand, their 
total catalog numbers would result in a different A^vc/A^WN ra- 
tio of 0.7, which might lead to biased estimates of the WR 
subclasses. iHamann et al.l (l2006l) selected only putatev ely single 
stars from the WR star catalog (Ivan der Hucht! 120011) for their 
study of the ratio between different spectral subclasses. The bi- 
narity status of invidual objects may have been revised in the 
meantime, and thus affect the number ratios. 
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For the Arches and the Central Parsec cluster the situation is 
similar, depending on the study, their field of view or the avail- 
able stellar subsample, the numbers i n the different column s of 
Table[T]may differ from those listed in Ivan der Huchd (l2006h . 

With these uncertainties in mind, we compare the number ra- 
tios Nwr/Nq and A^wc/A^vN with the prediction s of population 
synthe s is models as presented in th e figures by iLeitherer et al.l 
(Il999l) . IVazquez & Leithered (|2005 [ starburst99). The models 
account for so l ar me tallicity and single star evolution. From 
ILeitherer et al.l (1 19991 (Fig. 39 to 42)), we find lower Umits for 
the age of the Quintuplet cluster of 3.3 to 3.6 Myr in the scenario 
of instantaneous burst star formation. This agrees very well with 
the above HRD and isochrones. Although we cannot exclude 
a past supernova in the Quintuplet, the determined age would 
also be conform with the predictions of Leitherer et al.l (Il999l 
(Fig. 43)) in terms of first supernovae taking place only at about 
3.5 Myr in the star burst scenarios. 

One issue with the Quintuplet stellar population is obvi- 
ous: All WN-type stars in the Quintuplet have been classi- 
fied as WNL, i.e . they still contain hydrogen in their atmo- 
sphere (iLiermann et al. 2009, 2010). WNE stars (hydrogen-free 
WN stars) are significantly missing, although th e number ratio 
^^^M^wR is compatible with the value found by Ivan der Huchti 
(1200 lb for the Milky Way. From the regulai- Milky Way WR pop- 
ulation we woul d expect a sl ightly higher number of WNE stars 
than WNL stars (iHamann et al. 2006). WR 102c, which hes out- 
sid e the LHO field of view in the Sickle nebula, was classified 
bv lFiger et al.l (Il999ah as W N6 type star and thus b elongs to an 
"early" subtype (WNE). But iBarniske et al.l (l2008l) found from 
comparison with model spectra that the stellar atmosphere still 
contains up to 20 % hydrogen, making it a "WNL" star in the 
evolutionary sense. Interestingly, WNE stars seem to be miss- 
ing in the other clusters in the Galactic center region as well 
(see Liermann et al. 2010, and Table[TJ. From the Geneva mod- 
els we know that hydrogen burning lifetimes and WNL lifetimes 
are increased when the effects of rotation are included, while 
WNE lifetimes are hardly af fected. This might skew the number 
ratios in favor of WNL stars. Liermann et al.l(l2010l ) noticed that 
evolutionary models without rotation represent the Quintuplet 
WNL stars better than those with rotation. This might imply that 
the initial rotational velocity of 300km/s~' in the mo dels is too 
much. Or, as pointed out by IVanbeveren et al.l (l2007l) . the mis- 
match between evolution tracks and observations for massive 
stars may be due to inappropriate assumptions of the mass loss 
in the RSG stage. 

So far we have only considered single star evolution, but 
what about binaries? The most promising indications for bi- 
narity seem to be present among the WR stars, so we will 
focus on those in the following. In the Milky Way about a 
quart er of all WN stars a nd half of all WC stars are in bina- 
ries (Ivan der Hucht 20011 see als o Table[T}. As noted by sev 



eral authors (e . g. 



vanderHuch"3 l2006. Crow ther et all I2OO6I 



ILiermann et al.l |2009|) . TiT-band spectra that show an increas 
ing continuum with wavelength (IR excess) can be interpreted 
as being colliding w ind binaries in which dust is formed 
dWilliams et al.|[l987h . The majority of WC stars in current cat- 
alogs and surveys are binary candidates on this "dust" argument 
alone. 

With no time series of spectral observations available, the 
LHO spectra allow only indirect conclusions about binarity to 
be drawn. Among the WC stars in the Quintuplet are nine stars 
classified as WCd, i. e. dust producing colliding wind binary can- 
didates. Five of those, the enigmatic Quintuplet pinwheels, can 
be considered confirmed binaries as they were partially or fully 



spatially resolved by interferometrvbv rTuthill et al.l(l2006l) . Two 
further WC stars, LHO 76 and 79, have been classified in the 
LHO catalog as potential binaries on the basis of a "dusty" spec- 
trum. So, the number of WC binaries might even be higher 

Among five WN stars in th e LHO catalog, we find only 
one binary candidate: LHO 110. ILiermann et al.l (|2010) found 
small absorption lines in the spectrum which might belong to 
an O-type companion. This is corroborated by its high mea- 
sur ed radio flux. One fu rther WC and one WN star are listed 
by Ivan der Huchti (1200 6*) as binary candidate. But for these, 
WR 102f (WC8-(-?) and WR 102i (WN9+?), the LHO specti-a 
look unsuspicious, which is why the classification as binary 
candidates was dropped. The number of unconfirmed binaries 
makes it diflicult to compare the binary fraction with other 
Galactic center clusters and the Milky Way (see TablefTji. But 
it seems that the binary fraction of the Quintuplet WN and WC 
stars resembles roughly the one for the Milky Way in general. 

IVanbeveren et al.l (1 19981) argue that taking binaries into ac- 
count is required to get correct N^r/Nq and A^vc/A^WN ratios. 
Their population synthesis suggests that, especially for the inner 
Milky Way, a A^wr /Nq ratio of the order of 0.2, a value similar to 
the one found for the Quintuplet, cannot be explained with a con- 
tinuous star formation scenario. It could rather be explained with 
an enhanced star formation or star burst. In addition, according 
to those models a star burst that includes binaries encounters 
the effect of rejuvenation, i. e. mass transfer which takes place 
in massive close binaries leads to a population of young O-type 
stars mimicking a much younger age. This is supposed to happen 
after about 4 Myr, and would manifest in a group of blue young 
stars above the main sequence turn-off of the main population. 
We don't see a prominent group of such O-type stars in the HRD 
of the Quintu plet, which elim i nates a cluster age of > 6 Myr 
according to Vanbever en et al.l (1 19981 Fig. 1 1 B). Given the er- 
ror bars in our emp irical T^ff and log^ (Fig .Bli, and comparing 
this to the spread in IVanbeveren et al] (Il998l Fig. 11 A), we can 
furthermore limit the cluster age to < 4 Myr. This is within the 
range we estimate from the main sequence of OB and WR stars 
and the isochrones. 

F rom the good agree ment with the overall numbers ra- 
tios dvan der Huchti 1200 lb and the a ge hmits inferred frorn 
sing le and/or binary star evolution (IVanbeveren et al.l Il998l 
Leit herer et al . 1999) it seems secure to conclude that the for- 
mation of the Quintuplet cluster most likely needed a star burst 
event about 3.5 + 0.5 Myr ago. 



4.3. Late-type stars in the Quintuplet? 



According to IVanbeveren et al.l (1 19981) almost one third of all 
Galactic open cluster and stellar aggregates host WR stars and 
RSGs simultaneously. In the majority of them the total numbers 
are rather small, of the order of one or two, and in all cases the 
number of WR stars exceeds those of RSGs. 

In terms of abundance of these special stars, Westerlund 1 
(Wd 1) is the most prominent Galactic open cluster known to si- 
multaneously contain evolved early-type stars, such as luminous 
blue variables (LBVs), yellow hypergia nts (YHGs), WR stars , 
and RSGs. Based on single star evolution. lCrowther et al.l(l2006l) 
concluded from an number ratio of A^Rs r.+YHfi/A^vR = 8:24 a 
cluster age of 4.5-5 Myr for Wd 1. Meng el & Tacco ni-Garma^ 
(2007) list four RSGs plus six YHGs'found by Clark etal] 
(2005) w hich changes t he number ratio slightly (see Table[T}. 
Vanbeve ren et al.l (119981) argue on basis of single star lifetimes 
that WRs and RSGs can be present simultaneously only in a 
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Table 2. Comparison of stellar parameters for LHO 007. 
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Fig. 7. Detail of Fig. |4] with focus on the late-type stars with 
log(L/Lo) > 3.4 that might be Quintuplet RSG candidates. 
Tracks show stellar evolution models from iMevne t & Maeder 
( l2003h with and witho ut rotation (soli d and dashed lines, respec- 
tively) and from .Bressan et al.l (Il993l solid line with circles for 
Mini = 6 and 7 Mq). 



very limited time span of 4 to 5 Myr Too early no RSG would 
be present and too late the WR stars would be already gone. 

We already noted that in the Geneva mode ls with rota- 
tion th e WNL lifetime is increased. In addition, ^ Hirschi et al.l 
(|2004|) found that in the rotating models the RSG phase could be 
reached even before helium ignites in the stars, much earlier than 
in non-rotating models, and that the RSGs become more lumi- 
nous. Similar effects were found for the WNL stars. Both WNL 
and RSGs can appear spectra-wise as evolved stars, while they 
are still core-hydrogen burning. This might twist the number ra- 
tios towards more WNL than WNE stars plus the simultaneous 
presence of RSGs even at young cluster ages. 

For the Quintuplet cluster, eight late-type stars in the regime 
of log (L/Lq) > 3.4 are found (see Fig.Q. Including stellar 
evolution tracks from Bressan et al. (1993), most stars below 
log (L/Lq) < 4 could also be regarded as intermediate-mass stars 
(tracks with Mini ~ 6 to 7 Mq). However, only five out of eight of 
these stars are from the corrected sample of late-type stars, i. e. 
potential cluster members, due to their RV. This would result in 
a number ratio A^rsg+yhg/A^wr of 5:21 similar as for Wd 1 (see 
TablelTJ and could imply that the Quintuplet cluster is slightly 
older than just 4 Myr 

The stars LHO 006, LHO 036, LHO 053, and LHO 063 have 
not yet been studied in detail. The remaining stars are discussed 
individually in the following: 

LHO 007 - alias GMM7 (Q7), qF 192 

iGlass et aP d 19901) mentioned this promine nt and bright source 
as a late-type star in the Quintuplet region. iMoneti et al.l (1 19941) 
found the star to be an M-type supergiant with K - 7.36 mag, 
and Glass et al. ( 1999b) consider the star to be not variab l e. The 
stellar parameters from a detailed study by iBlum et all (l2003h 
are compared with our results in Table|2] Similar /T-band mag- 
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nitudes but slightly different effective temperatures put the star 
at the highest luminosity within our group of late-type stars. 
A number of studies also focused on the chemical composi- 
tion of the star (e . g. iRamirez et all l2003l ICunha et al.1 l2007l 
iDavies et al.ll2009l) . finding iron abundances to be nearly con- 
sistent with solar values but slightly enhanced a elements. Most 
authors agree that LHO 007 is an evolved star with an initial 
mass likely between 15 and 20 M© depending on the stellar evo- 
lution model. From our results we find that all stellar evolution 
models (with and without rotation, from Geneva and Padova) 
imply an initial mass of 15 M© (Fig.|7]l- The star's member- 
ship to the Quintuplet cluster was debated in terms of extinc- 
tion is sues (.Glass et al.1 1 19901 iMoneti et all 11994 iGlass et al.1 
Il999bl) . According to its T^V = 124kms"Mt may belong to the 
Quintuplet cluster. 



LHO 1 08 - ahas [NWS90] A, [GMC99] D 322, qF269 

This s tar is consid ered to be a foreground star by Nagat aet al.l 
(119901) . In contrast. [Glass et all (Ii999b.) fist the star in the group 
of probable cluster m embers. Its j?V = 127kms~' puts it in 
our corrected sample. ' Figer et all (Il999ab list the star in their 
Quintuplet sample, with a spectral classification as OBI and a 
luminosity of log(L/L0) = 5.54. But the LHO spectrum defi- 
nitely shows CO bands classifying it as late-type star (see LHO 
catalog spectra). Thus we derive a much smaller luminosity of 
log {L/Lq) = 3.5 1 . It places the star at the lower end of the group 
in vicinity to the non-rotating 9 M© track from Geneva or the 6 
to 7 Mo tracks from Padova (see Fig.|7). 



LHO 1 15 - alias GMM5 (Q5), MGM5-5, qF270N 
This star is one of the five from which the Quintuple t cluster de- 
rived its name. Its /T-band magnitude is 8.71 mag dGlass et al.l 
199 0). Later work s give a range o f values: K = 9 .1 mag 
(M onetietal.lll992h j i: = 8 78 mag (lMonetietalJll994. K = 
7.81 ma g (iFiger et al.iri999aL spectral type "late"'). iGlass et al.l 
(Il999bh list the star in their sample of long-period variables in 
the Quintuplet with an average K - 8.62 mag and a 680 d pho- 
tometric variability of 1.5 mag. They attribute the variability to 
the star being a Mira variable and argue on basis of stellar evolu- 
tion that the star would then be too old to be a Quintuple t mem- 
ber However, according to its photometric colors (Mone tiet al.l 
1994) and its apparent magnitude b eing in t he range expected for 
Miras in the GC region (.Glass et al]ll999ah . the cluster member- 
ship might still be likely. 

From our K - 8.6 mag we derive a luminosity of 
log(L/L0) = 4.38, which makes the star the second most lumi- 
nous in this group of RSG candidates. Its position in the HRD 
(Fig. IT) puts it at the stellar evolution track for a star with an ini- 
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Fig. 8. Luminosity function of the Quintuplet cluster for all stars 
in the LHO catalog (dashed line) and for the corrected sample 
(solid line). The photometric completeness at /^s - 13 mag is 
indicated by the vertical dashed line. The linear fits to both dis- 
tributions are indicated by the slope coefficient a. 



Fig. 9. IMF for the early-type OB stars in the corrected sample of 
the LHO catalog plus the WN stars. The distribution shown has 
variable mass range bins with a constant number of stars per bin 
to avo id numerical biases according to Mai z Apellaniz & TJbedal 
(l2005h : the linear fit gives a slope coefficient as indicated in the 
plot. 



tial mass of 9Mq. With the RV - 137kms"' it dropped out of 
the corrected sample, but still is within the 3 cr interval of possi- 
ble cluster membership. 



LH0 156 - alias [GMC99] D 307 

iGIass et al] ([l999b) list this star with K - 10.83 mag in the sam- 
ple of probable Quintuplet members and show a slightly variable 
light curve. However, no ffi'm conclusion about a possible vari- 
ability is drawn. We find K - 10.4 mag in the LHO catalog. It's 
derived luminosity log(L/L0) = 3.72 puts it on the non-rotating 
9 Mq track from Geneva or the tracks from Padua with 6 to 7 M© 
initial mass. However, with RY - 169 kms' the star is outside 
the limits for the corrected sample and even outside of the 3 cr 
cluster RY . 

Summarizing, there are a few KM stars in the Quintuplet field 
which share the radial velocity of the cluster One of them, 
LHO 007 alias Q7, has supergiant luminosity. Other ones are 
red-giant branch stars, assuming Galactic center distance. One 
can exclude that they are foreground KM dwarfs, because of 
their radial velocity and because of the '^CO signature in their 
spectra that is indicative for giants and supergiants. It would be 
very puzzling if these relatively old (>30Myr) stars would be 
gravitationally bound to the Quintuplet cluster Otherwise, they 
may belong to an old field population in the Galactic center re- 
gion. 

It seems that the assessment o f the cluster members hip, 
e. g. from proper motion studies (i HuBmann et al.l 1201 ll) . is 
needed to settle the question if older stars belong to the 
Quin tuplet. As discussed in the previous section, binary evolu- 
tion dVanb everen et al.lll998h and stellar evolution with rotation 
dHirschi e t al. 2003) may explain the coexistence of RSG and 
WR stars in a cluster of about 4 Myr, as the likely age of the 
Quintuplet. Alternatively, if none of these scenarios works, the 
conclusion would become inevitable that the cluster contains a 
second, older population of stars. 



4.4. Luminosity function and mass function 

To further exploit the LHO data of the Quintuplet cluster, we es- 
tablish the luminosity function by counting the number of stars 
per magnitude bin in our total cluster field of view; in addi- 
tion we apply a normalization factor to account for the cluster 
area. The resulting distribution is shown in Fig.|8] A power law 
with an exponent a = 0.24 + 0.06 can be fitted to the distribu- 
tion of the corrected sample. iFiger et all (Il999bl) argue that such 
kind of shallow distr i bution is expected for young coeval clus- 
ters. iHaravama et al.l (l2008l) found for example a - 0.27 for the 
young massive star-forming region NGC 3606. 

For a young coeval cluster it is expected that the mass dis- 
tribution in the stellar population still resembles the initial mass 
distribution. However, the evolved WR stars (WN and WC sub- 
types) and the presence of RSG candidates clearly show that 
some stars have already undergone significant mass-loss. On the 
other hand, we found that the majority of early-type stars form 
a main sequence around a 4 Myr isochrone. Thus, we assume 
the OB stars to have formed during one star formation event and 
omit the evolved WR stars and late-type stars for this moment 
from the further discussion. 

From the Geneva 4 Myr isochrone a mass-luminosity re- 
lation is constructed and applied. Our completeness limit of 
A^s - 13 mag translates to a minimum mass of Mini = 10 Mq 
for an OB star Thus we obtain initial masses in the range of 
10 to 78 Mq, with most stars concentrating in the 10 to 30 M© 
range and much fewer stars for higher mass bins. The error bars 
are estimated on the basis of the uncertainty of the luminosity 
propagating through to the derived masses. 

To determine the IMF we follow Mai z Apellaniz & Ubedal 
(l2005h and use bins with a variable mass range but a constant 
number of stars per bin to be robust against numerical biases. 
In Fig.|9]we show the resulting distribution for the OB stars in 
the corrected sample plus the WN stars with their masses as 
determined by Liermann et al. (2011). A linear fit of the form 
log(dA^/dM) ~ y X logM gives a power law with y = -1.23 + 
0.51. The fit for the OB stars only gives a slightly steeper result 
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of 7 = -1.66 + 0.51 as the WN stars basically populate the very 
high-mass bins. Both scenarios suggest a slightly top-heavy IMF 
in comparison to the canonical Salpeter-IMF (Salpeter 1955, y 
= -2.35). 

Tests with the Padova 4 Myr isochrone give similar results. 

Top-heavy IMFs are discussed for massive stellar cluster 
throughout the Gal axy and espec i ally fo r the Galactic center re- 
gion. For example, iManess et al.l (l2007h find j - -0.85 for the 
Central Parsec cluster assuming continuous star formation, albeit 
iPaumard et al.l (12006') favor a star burst scenario for this cluster 
The Arches cluster was analyzed by Stolte et al. (2005) who dis- 
cuss a flat present-day mass function with y = -1.9 to -2.1, but 
iPortegies Zwart et al.l (l2002i) do not exclude a regular Salpeter- 
IMF for that cluster when mass segregation is taken into account. 
Our results on the Quintuplet cluster seem to be in agreement 
with the findings of a recent study by iHuBmann et a l.i (.2 01 L) . 



5. Summary 

In this paper we have used the A'-band spectra of the LHO cata- 
log to obtain stellar temperatures and luminosities for all early- 
and late-type stars in the Quintuplet cluster Furthermore, we 
used the LHO catalog to obtain number ratios for the different 
stellar subclasses. As result we present a cluster HRD that is 
compared with stellar evolution tracks and isochrones; the main 
findings are the following: 

1 . The early-type OB stars in the Quintuplet are high-mass stars 
with initial masses > 8 M^. They form a main sequence 
population of massive stars around an isochrone of about 
4 million years. 

2. For the WN stars we determine an age of 3 million years 
fro m the isochrone s whic h is consistent with the ages found 
by iLiermann et al.l (l2010h . A possible interpretation for the 
slight diiference to the age of the OB stars is that the most 
massive stars might have formed late in the cluster formation 
process. 

3. Alternatively, if all WR stars in the Quintuplet evolved as 
secondary stars in binary systems, they might have been re- 
juvenated by mass transfer while being originally older 

4. Considering different number ratios (e.g. N^/r/Nq, 
A^vc/A^vN) in the Quintuplet we find best agreement from 
comparison with population synthesis models for a star burst 
that formed the Quintuplet cluster about 3.5+0.5 million 
years ago. 

5. The late-type KM stars need much longer to evolve to their 
present stage. According to their radial velocities the major- 
ity of them cannot be considered as cluster members. Five 
of the KM stars, however, share the radial velocity of the 
Quintuplet cluster One of them (LHO 007 alias Q 7) is a su- 
pergiant, the others are less luminous. They may represent 
an older field population in the Galactic center region. 

6. If LHO 007 really is a physical member of the Quintuplet 
cluster, this would have interesting implications. With stan- 
dard single-star evolution, the simultaneous presence of WR 
stars and a RSG of initially about 15 Mq is only expected for 
clusters of much older age. 

7. Recent stellar evolution models which take the effects of ro- 
tation into account, as well as considering binary evolution 
channels, may help to explain the coexistence of WR stars 
and RSGs in a cluster of only about 4 million years. 

8. We derive the luminosity function of the Quintuplet cluster 
with a power-law slope of a = 0.24 + 0.06. This is consis- 



tent with a young coeval cluster and confirms the star burst 
scenario for the formation of the Quintuplet. 
9. For the subgroup of the high-mass stars we derive a mass- 
luminosity relation from the 4 million year isochrone and 
apply it to obtain the IMF. The result is a power law with 
a slope that suggest a slightly top-heavy IMF compared to 
the Salpeter-IMF If we include the WN stars the eff'ect is 
more pronounced. 

In summary, the Quintuplet cluster seems to have an interme- 
diate age of 3.5+0.5 million years compared to the other Galactic 



center clusters, the Arches (2.5+0.5 million years, Martins et aL 
2008) and the Central cluster (6 ±2 million y e ars, P aumard et aL 
2006). As recently discussed bv lBoker et al.l (l2008l) . an age gra- 
dient of young stellar clusters in extragalactic nuclei can be ex- 
plained in a scenario of almost continuous star formation by 
over-density regions in a circum-nuclearring (CNR) in a galaxy. 
This would lead to a burst-like release of stellar clusters in reg- 
ular time steps like "pearls on a string" in the central zone. 
From radio data we know the existence of the Milky Way's 
CNR and the active high -mass star formation region SgrB2 (e. g. 
Morris & Serabvn 1996). Thus, if the star formation in the inner 
Milky Way is interpreted in the way described above, this would 
explain the observed age sequence in the three Galactic center 
clusters. 
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